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We study the potential of International Linear Collider (ILC) at yfs = 500 GeV to probe new dark 
matter motivated scenario where the bottom squark (sbottom) is the next to lightest supersymmetric 
particle. For this scenario, which is virtually impossible for the LHC to test, the ILC has a potential 
to cover a large fraction of the parameter space. The challenge is due to a very low energy of jets, 
below 20-30 GeV, which pushes the jet clustering and flavour tagging algorithms to their limits. The 
process of sbottom pair production was studied within the SiD detector concept. We demonstrate 
that ILC offers a unique opportunity to test the SUSY parameter space motivated by the sbottom 
neutralino co-annihilation scenario in cases when the sbottom production is kinematically accessible. 
The study was done with the full SiD simulation and reconstruction chain including all Standard 
Model and beam backgrounds. 
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I. INTRODUCTION 

Among the best candidates for a theory beyond the 
Standard Model (SM), Supersymmetry (SUSY) [IHl] re- 
mains a very compelling theory even after 30 years with- 
out an experimental confirmation. SUSY is very attrac- 
tive because it successfully solves principal theoretical 
and experimental problems of the SM. Supersymmetric 
theories provide a natural solution to the gauge hierar- 
chy problem of the SM and incorporate the unification 
of gauge coupling constants. Furthermore, the lightest 
supersymmetric particle (LSP) is stable if the R-parity 
is conserved and can serve as a good cold dark matter 
(CDM) candidate. Besides, SUSY has all ingredients to 
provide a solution to the baryogenesis problem via the 
intermediate scale leptogenesis [7| or via the electroweak 
baryogenesis [8|, [9j. 

The very existence of CDM is a crucial argument in 
favour of SUSY, and, at the same time, an important 
SUSY constraint. The most direct evidence for CDM in 
the Universe comes from observations of galactic rotation 
curves. Binding of galaxies in clusters, matching obser- 
vations of large scale structure with simulations, gravita- 
tional microlensing, baryonic density of the Universe as 
determined by Big Bang nucleosynthesis, observations of 
supernovae in distant galaxies, as well as measurements 
of anisotropics in the cosmic microwave background radi- 



ation (CMB) can also be considered as strong confirma- 
tions of CDM (for reviews see e.g. [13, El)- I n particu- 
lar, the analysis of the Wilkinson Microwave Anisotropy 
Probe (WMAP) and galaxy survey data puts the most 
stringent constraint on the ratio of the dark matter den- 
sity to the critical density (l2j : 



VcDAih 2 = 0.11110.016 ( at 95%CL) 



(1) 



where h = 0.74±0.03 is the normalized Hubble constant. 

In most of the parameter space of SUSY models, the 
value of QcDMh 2 is well above the WMAP bound. In 
particular, in case of well explored minimal supergrav- 
ity (mSUGRA) model [13l - [la |. which is defined by uni- 
versal soft SUSY breaking scalar masses (mo), gaugino 
masses (m 1 / 2 ) and A-terms (Ao) at GUT scale, the CDM 
is the lightest neutralino Xi- Neutralino should annihi- 
late or co-annihilate with other SUSY particles at the 
Early Universe time intensively enough to lower down 
^CDAih 2 — il^ah 2 to experimentally acceptable level. 
This happens in the following special regions of SUSY 
parameter space: 

1. The bulk annihilation region at low values of mo 
and mi/2, where neutralino pair annihilation occurs 
at a large rate via i-channel slepton exchange. 

2. The stau co-annihilation region at low mo where 
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in the early universe [17|, \li 



3. The hyperbolic branch/focus point (HB/FP) re- 
gion [I9U22I ] at large mo near the boundary of 
the Radiative Electroweak Symmetry Breaking 
(REWSB) excluded region where the superpoten- 
tial Higgsino mass term becomes small and the 
neutralinos have a significant higgsino component, 
facilitating their annihilations into WW and ZZ 
pairs. 

4. The A-annihilation funnel, which occurs at very 
large values of tan j3 ~ 45 — 60 [23T - [28j . In this case, 
one has ttia ~ 2m^o . An exact equality in the 
mass relation is not necessary, since the A width 
can be quite large (Ta ~ 10 — 50GeV); one can 
thus achieve a large X1X1 ~> A ~~ ^ / '/ annihilation 
cross section even if 2m^o is several values of width 
away from resonance. The heavy scalar Higgs H 
also contributes to the annihilation cross section. 

In addition, for a particularly large Aq values there 
exists a region of neutralino top-squark (stop) co- 
annihilation and, at low TO1/2 values, a light Higgs h 
annihilation funnel region. 

In all these regions the processes of neutralino 
(co-) annihilation have a high enough rate to suppress the 
value of f^o h 2 to the experimentally acceptable level. 

One should note that each particular DM motivated 
region of SUSY parameter space defines a specific col- 
lider phenomenology. It is also worth to mention that 
some of those regions are problematic for the LHC and 
the ILC would play a crucial role in exploring them or 
even in discovering Supersymmetry in these regions. The 
representative example of such a region is the HB /FP re- 
gion which, as it was shown, can be efficiently covered 
only by the 1 TeV ILC collider [H-dTJ . 

In this paper we study the sbottom-neutralino co- 
annihilation (SBC) scenario and its ILC phenomenol- 
ogy. In the SBC scenario, as follows below, the typi- 
cal sbottom-neutralino relative mass difference <5rog~o = 

(rrc.£ — to^o )/m^o is below 10%. This parameter space 
is inaccessible neither at Tevatron nor at the LHC. It is 
interesting to notice that Tevatron would be in a better 
position than the LHC to study light sbottom scenario 
which requires very low px 6-jet threshold and would be 
sensitive to kinematically allowed parameter space with 
5mi^o ^ 40% [32j,l33j. Therefore the SBC scenario which 
we study here is new for the collider phenomenology. 
Also, as we show below, the SBC scenario has not been 
studied in detail from the theoretical point of view. 

The rest of the paper is organized as follows. In Section 
II we discuss in detail the SBC scenario and, in particular, 
the GUT scale boundary conditions providing SBC at the 
Electroweak Scale. In Section III we present the signal 
and background studies at the detector level and work 
out a strategy for observation of the signal. Conclusions 
are drawn in Section IV. 



As we mentioned above, one of the important mecha- 
nisms of neutralino relic density suppression occurs when 
the next-to- lightest supersymmetric particle (NLSP), 
is close in mass to the neutralino. In this case, 
the neutralino relic density is not only suppressed by 
the neutralino-neutralino annihilation, but also by co- 
annihilation with the NLSP. 

In the traditionally well-explored mSUGRA model 
only stau or stop co-anninhilation can take place to sup- 
press the CDM relic density down to the experimentally 
allowed region. 

On the other hand in the various models motivated 
by SUSY GUTs one can expect different kinds of devi- 
ation from universal boundary conditions which would 
predict qualitatively different phenomenology still con 
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sistent with the present experimental data (see e.g. 
and references therein). 

In such models some new co-annihilation scenarios 
can take place. In particular, in (35l . |36| it was no- 
ticed that under certain relation of scalar masses the 
sbottom-neutralino relative mass difference Smy^o = 

(m.£ — m^o)/m^o can become below of about 10% and 
provide effective SBC suppression of the relic density 
down to the experimentally acceptable level. The pa- 
pers (35I |36| concluded that in case of minimal sfermion 
non-universality (mSFNU) defined by (raio, TO5, Aq, 
TO1/2, tan /?, sign(fj,)) parameter space the SBC sce- 
nario could be realised, where the left-right squark and 
sfermion mass parameters at the GUT scale extend the 
universal mo mSUGRA parameter. 
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no — m Q = " l u — " l E 

It was found that the ratio 

K = m 5 /m w 

below of about 0.5 could provide a SBC scenario. 

In this section we take a closer look at the results 
of [H, Hi| where conditions for SBC were derived us- 
ing simplistic one-loop Renormalisation Group Equations 
(RGEs). In our paper we use mass spectra calculation 
from ISAJET v7.79 package [37} where the complete two- 
loop RGEs for the gauge couplings, Yukawa couplings 
and soft breaking terms are implemented. Another rea- 
son for checking of the SBC scenario was app lying a dif- 
ferent approach in comparison with (35l . l36l |. where the 
exact b — r unification was imposed at the GUT scale and 
RGEs were run just once from this scale. In the ISAJET 
the RGEs are solved iteratively, running from the weak 
scale to the high scale and back, while Yukawa couplings 
are not fixed at the GUT scale. After each iteration the 
SUSY masses are recalculated, and the renormalization 



3 



group improved one-loop corrected Higgs potential is cal- 
culated and minimized. This way of mass calculation is 
widely accepted to be the most accurate one. 

We have performed an exhaustive scan of the SUSY 
parameter space defined by 

< m 1/2 < 2TeV, |A | < 3TeV, 5 < tan/? < 50 
< m 5 < 5 TeV 
< m Q < 5 TeV 
< mu = m E < 5 TeV (2) 

where we went beyond the minimal sfermion non- 
universality (mSFNU) scenario, defined in [35|,[36j], and 
split mio = niQ — mu — mE into two parameters tuq 
and mjj — the while keeping m§ = mt = mo as defined 
previously. 

Our first results are shown in Figure [TJ which presents 
the results of the scan in the SUSY parameter space de- 
fined by (J2J) shown in mo/my vs. mQ/mu (left) and mQ 
vs. mQ/mjj (right) planes respectively with the follow- 
ing color code for Sm^-a = (rag — m^o)/m^o: 

a) 1.0 < 5rag- (red/dark), 

b) 0.2 < 6m~ b ~o < 1.0 (pink/grey), 

c) 0.1 < Sm^o < 0.2 (light-green/very light grey), 

d) Sm^-o < 0.1 (dark-green/light grey). 

One can see that the SBC scenario cannot be re- 
alized in the mSFNU scenario adopting a more real- 
istic approach which we used in our study. As indi- 
cated in Figure [TJ the neutralino-sbottom mass split 
5mj-o < 0.2 which is required by the SBC scenario, de- 
fines a non-universal condition mQ < O.bmu which even- 
tually cannot be realized in the mSFNU scenario with 
mQ = mu. We have also found that in the mSFNU re- 
gion for the SBC scenario from [3^, H|| , actually, the stau 
co-annihilation takes place when the ISAJET approach is 
used for the mass spectrum calculation. We illustrate this 
in Figure [5] which presents sbottom and stau masses ver- 
sus K = m^/mio for the mSFNU benchmark points from 
the mass spectra in (3|| for which was calculated using 
ISAJET v7.79. One can see that in the nSFNU case the 
stau-neutralino mass split becomes small for sufficiently 
small K, providing a stau-coannihilation mechanism of 
CDM suppression while sbottom is quite far from being 
NLSP. 

We have also found that if one departs from the 
mSFNU scenario and consider mg as an additional in- 
dependent parameter as we did in our scan, there are 
solutions for the SBC scenario. For sufficiently small val- 
ues of mQ < 0.5mu as shown in Figure [3] the bottom 
squark becomes NLSP with a mass quite close to the 
neutralino mass which provides the sbottom-neutralino 
co-annihilation suppression of CDM. Hereafter we refer 
to this scenario as next-to-minimal sfermion mass non- 
universality (nmSFNU). 

Finally, in Figure [T] (top) we present a mass spec- 
tra for the nmSFNU scenario for the benchmark point 



(toi/2 = 0.5 TeV, tan/3 = 38, A = ITcV and m v = 
m L = 3.95 TeV) relevant to the 500 GeV ILC which 
will be able to produce sbottom with mass below 250 
GeV with about 10% relative mass split with neutralino. 
For a sufficiently small toq parameter (of the order of 
1 TeV) the SBC scenario takes place. In Figure [4] (bot- 
tom) we present the respective behaviour of ilh 2 versus 
K = m^/mjj which reaches WMAP constraints when 

Sm bx° - 0A - 

We would like to stress that deviations from mSMNU 
which, as we have discussed, provide a sbottom co- 
annihilation scenario would be an indicator of non-trivial 
D-term contributions to soft scalar masses, which in its 
turn could give a hint about the gauge group content 
at very high energies as noticed in [381 ] . As we will show 
below, the ILC has a promising potential to cover param- 
eter space motivated by the SBC scenario thus shedding 
light on the gauge group properties at the GUT scale. 

III. SBOTTOM PRODUCTION AT THE ILC 
A. The ILC and the SiD Detector Concept 

The ILC is a future electron-positron collider designed 
to collide particles at the centre of mass energy of 
500 GeV. At the ILC sbottom pairs can be produced via 
s-channel photon or Z-boson exchange followed by sbot- 
tom decay into x? and b-quark, see Figure [5] In the sce- 
nario described above such events will result in two soft 
b-jets, which are generally not acoplanar nor acolinear, 
and a missing energy. 

The sensitivity to the sbottom production was studied 
in the framework of the Silicon Detector (SiD) concept 
using full detector simulation and event reconstruction. 
SiD is a detector concept [39[ designed for precision mea- 
surements of a wide range of possible new phenomena at 
the ILC. It is based on a silicon pixel vertex detector, sil- 
icon tracking, silicon-tungsten electromagnetic calorime- 
try, and a highly segmented hadronic calorimetry. Parti- 
cle Flow Algorithm (PFA) approach [39] is an important 
strategy driving the basic philosophy and layout of the 
detector. SiD also incorporates a 5T solenoid, iron flux 
return and a muon identification system. A schematic 
view of SiD quadrant is shown in Figure [6] 

B. Signal and Background Rates and Properties 

In this study the CalcHEP package f4C| was em- 
ployed to generate SUSY signal events. The events were 
converted to the Les Houches format [4l[, passed to 
Pythia [42j for their fragmentation and particle decays 
and consequentl y to the SiD full simulation and recon- 
struction chain [39j. Several points in the MSSM pa- 
rameter space were chosen close to the kinematic limit in 
order to investigate the discovery potential. These points 
correspond to various masses of the bottom squark and 
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FIG. 1: The results of the scan in SUSY parameter space defined by @ shown in [mo / mu-mQ /mjj) (left) and (mQ-rriQ/mu) 
(right) planes respectively with the following color code for Srrifca = (m^ — m^o)/m^o: 

a) 1.0 < Sm^o (red/dark), 

b) 0.2 < Sm-^o < 1.0 (pink/grey), 

c) 0.1 < Sm^o < 0.2 (light-green/very light grey), 

d) 5m- b ^a < 0.1 (dark-green/light grey). 
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FIG. 2: Sbottom and stau masses versus K = ms/mio for 
mSFNU benchmark points from [36l ] mass spectrum which 
was calculated using ISAJET v7.79. 



FIG. 3: Sbottom and stau masses versus K = ms/mu for 
mSFNU extended by m,Q as an additional independent pa- 
rameter which provides the SBC scenario 



neutralino: (230, 210); (240, 210); (230,220); (240,220); 
(220,210), all in GeV. For each signal point 200k events 
were generated, also accounting for the ISR, FSR and 
Beamstrahlung (BS) processes. 

The cross section of sbottom pair production depends 



strongly on how close the sbottom mass is to the ILC 
kinematical limit, see Figure [7] for y/s = 500 GeV, diving 
quickly down for masses approaching 250 GeV. The mass 
of the lightest bottom squark mass, m^, and the mix- 
ing angle of the lightest and heaviest bottom squarks, 9b, 
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FIG. 4: Sample point for nmSFNU scenario which is relevant 
for 500 GeV ILC: m 1/2 = 0.5 TeV, tan/3 = 38, A = 1 TeV 
and mu = tul ~ 3.95 TeV. For sufficiently small tuq param- 
eter (of the order of 1 TeV) the SBC scenario takes place. 
Top: sbottom mass versus K = m^/mu; Bottom: fi/i 2 versus 
K = ms/mu which satisfies the WMAP constraints when 



K « 0.555 and 8m 
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FIG. 5: Feynman diagrams for sbottom pair production (left) 
followed by sbottom decay into x? an d b-quark (right) 



completely determine the cross section of pair production 
of lightest bottom squarks. For the fixed mass, the 
cross section can vary by factor of about 5, depending on 
the value of 0b . The solid line in Figure [7] corresponds to 
the case of the minimal cross section, which takes place 
at cos0 7 b nm ~ 0.35, while the dashed line corresponds to 
the case of the maximal cross section which takes place at 
cos0 b nax = 1. The line grey levels (light/dark) in Figure[7J 
present the effect of Beamstrahlung radiation (on/off), 
which decreases the cross section by about 20%. The 
signal cross section is about 4 — lfb for 215 — 230 GeV 
bottom squark mass for 0™ Mn (solid dark line, Figure UJ) 
and about factor of 5 larger for Q™ ax (clashed dark line, 
Figure [JJ . In our analysis we use signal cross section 
0mm ^ therefore conservatively estimating the ILC poten- 
tial to probe the SBC scenario. 

While the signal is at a few fb level, the main back- 
ground, which originates essentially from quasi-real pho- 
tons emitted by electron and positron at very low scat- 
tering angles, 77 — > bb, is about 4.5 pb at 500 GeV 
ILC. Another significant background of similar topolog- 
ical properties is the di-jet background from e + e~ — > bb 
which has about 2.7 pb production cross section as well 
as e + e~ — > qq, where q represents other quarks misiden- 
tified as b. One of the powerful variables, which could 
suppress such dominant backgrounds, is the b-jet sepa- 
ration variable AR^b = 



A 2 



where A,* and A 



v 



is the difference of jet azimuthal angles and rapidities, 
respectively. In the case of backgrounds with two back- 
to-back b-jets ARbb > 7T while in the case of signal, which 
has also two neutralinos in the final state together with 
two b-jets, ARbb can be well below tt. This is illustrated 
in Figure [S] 

Though the parton level ARbb distribution for signal 
versus background looks very promising and suggests to 
eliminate the e + e~ — > bb and 77 — > bb backgrounds com- 
pletely by setting ARbb > ^ cut, the situation is not that 
optimistic at the full detector simulation level: ARbb dis- 
tribution from 77 — > bb and e + e _ — > bb processes can 
"leak" into ARbb < n region due to detector resolution 
and event reconstruction effects. The realistic ARbb dis- 
tribution after the full detector simulation is shown as 
filled histograms in Figure |U 

Complete sets of unbiased Standard Model processes at 
500 GeV were generated using the Whizard Monte Carlo 
program [43| with a total statistics of about 7M events. 
All 0, 2, 4 and 6 fermion final states were generated, 
see |1 H for details. PYTHIA g2] was used for final 
state QED and QCD parton showering, fragmentation 
and decay to provide final-state observable particles. In- 
cluded in this sample are backgrounds arising from in- 
teractions between virtual and beamstrahlung photons. 
Both signal and background samples are considered un- 
polarised. The detector simulation is based on Geant 4 
toolkit [45|, Hlj]. A thin layer of Linear Collider spe- 
cific code, SLIC [13], provides access to the Monte Carlo 
events, the detector geometry and the output of the de- 
tector hits. 



6 



7000 



6000 



5000 



4000 



5000 — 



2000 — 



1000 




QDO 



Q 





* ' * » I ' ■ ij>' i > i t i ■ n ■ I i ■ i < i ■ r t « 



' 



1000 2000 3000 4000 5000 

FIG. 6: A plane view of a quadrant of SiD. Dimensions are in mm 



> l I I ■ I t I I I I l ■ I I ■ I I I L t l 

6000 7O0Q 8000 9000 



C. Analysis and Results 



The vertexing package used for the jet flavour tag- 
ging in this analysis was developed by the LCFI collab- 
oration [48[ . The main algorithm of the package, the 
topological vertex finder ZVTOP, reconstructs vertices 
in arbitrary multi-prong topologies. It combines relevant 
variables in neural networks which are then separately 
trained on samples of b-, c- and light quarks. The most 
discriminating variables are the corrected vertex mass, 
the joint probability, the secondary vertex probability, 
the impact parameter significance of the most significant 
track and number of vertices in the event. The joint prob- 
ability is defined as the probability for all tracks in a jet 
to be compatible with hypothesis that they originate at 
the primary vertex. The performance of the LCFI pack- 
age optimised for the SiD detector is shown in Figure [3] 
which shows the flavour tagging purity as a function of 
efficiency for the provided three tags. 



The SiD specific optimisation was performed by build- 
ing new neural networks and training them using di-jet 
samples which passed through the full SiD simulation and 
reconstruction. 

Observation of the sbottom production is challenging 
due to a very low energy of b-jets, below 20-30 GeV, 
which pushes the jet clustering and tagging algorithms 
to their limits. Figure [TU] illustrates that while the b- 
tagging efficiency is about 75% above 60 GeV, it is falling 
steeply for lower jet energies. The energy of b-quarks is 
determined by the mass difference between m~ h and m^o . 

Jets were reconstructed by employing the Durham kr 
algorithm with fc™ m = 10 GeV and limiting the max- 
imum number of jets to two. Figure [TT] shows distri- 
butions of several variables used in the analysis for the 
signal with m~ b — 230 GeV and m^o = 210 GeV and com- 
pares them to the inclusive SM background. The signal 
was scaled up by a factor of 10 5 for the presentation 
purposes. The variables are (from top-left): distance be- 
tween jets AR in r\ — cf> plane, maximum of jet rapidities 
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FIG. 7: Sbottom pair production cross section at ILC at 
y/s = 500 GeV: solid line corresponds to the case of the min- 
imal cross section which takes place at cos8™ m ~ 0.35 while 
the dashed line corresponds to the case of the maximal cross 
section which takes place at cos9™ ax — 1. The line (light/dark 
grey) colors present the effect of Beamstrahlung radiation 
(on/off). 
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FIG. 8: 
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e+e~ -> bibi -> fefcXiX? signal (230,210) (empty histo gram 
with dotted line border) and e + e~ — > bb background (empty 
histogram with dashed line border) processes compared to 
the full detector simulation level for signal (light grey shaded 
histogram) and background (light grey shaded histogram). 
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FIG. 9: Flavour tagging performance of the LCFI package in 
terms of purity vs. efficiency optimised for the SiD detector 
concept and evaluated for di-jet events e + e~ — > qq. Three 
tags provided are shown: b-tag, c-tag and c-tag on the b- 
background only (c(b)-tag). 
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FIG. 10: Jet b-tagging efficiency as a function of the jet energy 
for b-jets from di-jet events e + e~ — ► qq. 



mace ( 1 771 1 , 1 772 1 ) , visible energy, acoplanarity, number of 
particles in event, number of charged particles, number 
of identified leptons, number of particles in the forward 
calorimeter used for veto, neural nets outputs of jets for 
b-tag, c-tag, c(b)-tag and the momentum isotropy. 

Majority of the SM background, for the exception 
of the two-photon background, is suppressed by the 
E V isibie < 80 GeV selection. The signal events are mis- 
balanced in energy due to the neutralinos. The dom- 
inant two-photon and e + e~ — ► qq backgrounds pro- 
duce jets in a back-to-back topology motivating the 
AR v ,p < 3.0 selection. To further suppress the back- 



8 



ground max(\r]i\, |f^|) < 2.0 is required, where rji and 
r]2 are jet pseudorapidities. Finally, the total number 
of reconstructed particles in an event is required to be 
within 10 < N par ticies < 60 for the mass point (230,210) 
and it is adjusted for the other mass points individually. 

An important ingredient in this analysis is an elec- 
tromagnetic veto using the forward detector at very low 
polar angles, above 9 = lOmrad. This is used to sup- 
press the two-photon background as well as other back- 
grounds with a large ISR contribution. The forward de- 
tector acceptance was estimated using a simple geomet- 
rical model. If a photon or electron with E > 300 MeV is 
detected within the acceptance the event veto is applied. 

For the final event selection a FANN neural net- 
work [49| based on the above variables was defined and 
trained on independent signal and background samples. 
The resulting neural net output is shown in the left plot 
of Figure [TJl The result is interpreted in terms of signal 
significance calculated as S/ y/ S + B, where S is the num- 
ber of selected signal events and B the number of selected 
background events. The event numbers are normalised 
to the total luminosity of 1000 fb _1 . Figure H"2l shows a 
distribution of S/y/ S + B as a function of the number 
of selected signal events with each bin corresponding to 
a particular selection of the final cut on the neural net 
classifier output. Big variations are caused by the back- 
ground events with large weights. The large weights were 
used for some SM processes with large cross sections. 

Based on the above the signal cross section statisti- 
cal uncertainty was calculated to be equal to 15 % for 
the (230,210) mass point. The mass points (240, 210); 
(220,210); (240,220) all can be observed with 3 a statis- 
tical significance and the mass point (230,220) with 2 a 
statistical significance. A conservative estimate of the 3 a 
observation region is shown in Figure 1131 



Such a measurement, however, relies on a presence of 
the forward detector used to veto two-photon background 
events and soft di-jet events with a dominant ISR contri- 
bution. 

The approximate reach of the 500 GeV ILC extends to 
within 10 GeV to the kinematic limit for the sbottom 
pair production and within 5 GeV to the kinematic limit 
of the sbottom decay to LSP and b-quark. 
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IV. SUMMARY 



We have performed analysis of ILC potential to probe a 
sbottom co-annihilation (SBC) scenario. We have found 
that if the SBC scenario is confirmed at the ILC, this 
would be an indication of related non-minimal sfcrmion 
non-universality scenario at the GUT scale reflecting 
non-trivial D-term contributions to soft scalar masses. 
This would give a hint about the gauge group content at 
very high energies. 

The probe of SBC scenario is very challenging even at 
the ILC due to the fact that background is higher than 
the signal by several orders of magnitude. One should 
stress, that the realistic detector simulation and recon- 
struction is crucial for this study due to the important 
role of resolution effects. 

Using neural network analysis with the optimised in- 
put variables we conclude that the SBC scenario can be 
probed very close to the ILC kinematic limits and, in the 
case of discovery, the sbottom production cross section 
can be measured with a reasonable statistical precision. 
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FIG. 11: Neural net input variables. Signal with = 230 GeV and m^o = 210 GeV (open histogram) is shown together with 
the inclusive SM background (filled histogram). The signal was scaled up by a factor of 100000. The dashed lines indicate 
selection cuts applied before neural net training, see text for details. 
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FIG. 12: Neural network output (left) is shown for events which passed the basic selections for signal with = 230 GeV and 
m^o = 210 GeV (open histogram) and inclusive SM background (filled histogram). Statistical significance (right) is shown as a 
function of a number of selected signal events. Dashed line shows a limit of four standard deviations. Big variations are caused 
by SM events with large weights. 




FIG. 13: Selected mass points in the - m^,o plane. The 
dashed line corresponds to the kinematic limit corresponding 
to a small mass difference between my and m^o. The filled 
area is an estimate of the 3<r observation region. 
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